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Abstract 

We report spin wave and DMRG studies of the ground and low-lying excited 
states of uniform and dimerized alternating spin chains. The DMRG procedure is 
also employed to obtain low-temperature thermodynamic properties of the system. 
The ground state of a 2N spin system with spin-1 and spin-^ alternating from site 
to site and interacting via an antiferromagnetic exchange is found to be ferrimag- 
netic with total spin sq = N/2 from both DMRG and spin wave analysis. Both the 
studies also show that there is a gapless excitation to a state with spin sq — 1 and 
a gapped excitation to a state with spin sq + 1- Surprisingly, the correlation length 
in the ground state is found to be very small from both the studies for this gapless 
system. For this very reason, we show that the ground state can be described by 
a variational ansatz of the product type. DMRG analysis shows that the chain is 
susceptible to a conditional spin-Peierls' instability. The DMRG studies of mag- 
netization, magnetic susceptibility (%) and specific heat show strong magnetic-field 
dependence. The product \F shows a minimum as a function of temperature (T) at 
low-magnetic fields and the minimum vanishes at high- magnetic fields. This low- 
field behaviour is in agreement with earlier experimental observations. The specific 
heat shows a maximum as a function of temperature and the height of the maxi- 
mum increases sharply at high-magnetic fields. It is hoped that these studies will 
motivate experimental studies at high-magnetic fields. 



1 Introduction 

The low-dimensional magnetic systems exhibit a wide variety of exotic physical phenom- 
ena. Many of these interesting and novel phenomena were first predicted from theoretical 
studies on spin systems in one-dimension [I], 0, [|. Fascinating amongst these have been 
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the spin-Peierls instability and the Haldane conjecture. These predictions motivated a 
number of experimental efforts towards synthesis of low-dimensional magnetic systems. 
The experimental measurements @, [Tj on some of these systems have provided support 
for the Haldane conjecture and also for the existence of spin-Peierls instability in quasi- 
one-dimensional systems. While experimental studies^, || of the spin-1 antiferromagnet 
Ni(C2H%N2)2N02{Cl04) clearly show the existence of Haldane gap, the measurement of 
the magnetic properties in a series of quasi-one-dimensional compounds confirm the pres- 
ence of spin-Peirels instability at low temperatures 0. In recent years, purely inorganic 
systems that show low-dimensional behaviour have also been synthesized. The compound 
CuGeOj, has been shown to be a quasi-one-dimensional system exhibiting spin-Peierls 
instability 0. Another class of compounds that have become important in recent years 
are the systems with spin ladders. The compound that closely approximates a spin ladder 
is vanadyl pyrophosphate with molecular formula (VO)2-P207- A spin ladder with hole 
doping is predicted to support hole binding based on exchange energy considerations as 
two holes occupying the same rung of the ladder would be energetically favoured, if the 
exchange constant for the rung is larger than that along the leg. These systems are hence 
expected to exhibit non-BCS type of superconductivity |§. 

Yet another challenge in the area of molecular magnetism has been to synthesize a 
molecular ferromagnet. While there have been many models for ferromagnetic exchange 
in molecular systems ||, the actual synthesis of molecular ferromagnets has been rela- 
tively recent. The search for molecular ferromagnet has led to the discovery of many 
interesting molecular magnetic systems. In recent years, quasi-one-dimensional bimetal- 
lic molecular magnets, with each unit cell containing two spins of different spin value 
have been synthesized [ fTD] . These systems contain two transition metal ions per unit cell 
and have the general formula ACu(pbaOH)(H20)^.2H20 with pbaOH=2-hydroxo-l,3- 
propylenebis(oxamato) and A = Mn, Fe, Co, Ni and belong to the alternating or mixed 
spin chain family These alternating spin compounds are seen to exhibit ferrimagnetic 
behaviour. The ferromagnetic or ferrimagnetic alignment of spins in molecular systems 
is usually difficult to achieve due to the diffused nature of the molecular orbitals (MOs) 
which the unpaired electrons occupy. The direct exchange integrals involving MOs are 
much smaller than that found in transition metal ions. However, even with the smaller 
exchange integral, we could have a ferro or ferri magnetic spin alignment if the molecules 
are so arranged as to yield a small intermolecular transfer integral. In such a situation, the 
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kinetic stabilization of the electrons is rather small and parallel alignment of the effective 
spin of the electrons on the molecules is favoured leading to a magnetic ground state |T^| 

Thermodynamic properties of these alternating spin compounds show very interest- 
ing behaviour [|TT|, [13| . In very low magnetic fields, these systems show one-dimensional 
ferrimagnetic behaviour. The %T vs T (x is the magnetic susceptibility and T is the tem- 
perature) plots show a rounded minimum. As the temperature is increased, xT decreases 
sharply, goes through a minimum before increasing gradually. The temperature at which 
the minimum occurs differs from compound to compound. The behaviour of field induced 
magnetization with temperature is also quite interesting as the ground state of the system 
is a magnetic state. At moderate magnetic fields, with increase in temperature, the mag- 
netization slowly increases, shows a broad peak and then decreases. Such behaviour is 
conjectured to be due to irregular variation in the spin multiplicities of the energy levels of 
the system. Theoretically, there has been an earlier study of the alternating spin sytem in 
three dimensions that focuses on the dependence of Curie temperature on the anisotropy 
of the spins [|l4l| . 

Motivation for the present study comes from the above experimental observations on 
the quasi-one-dimensional alternating spin systems. Theoretical studies of spin chains 
so far, have however been concerned mainly with the antiferromagnetic spin systems 
with unique site spin, exceptions being dilute spin impurity problems. In this study, 
we first analyse the ground and low-lying excited states of the mixed or alternating an- 
tiferromagnetic spin chain by employing a spin wave theory. The spin wave theory is, 
however, not as accurate as the recently developed density matrix renormalization group 
(DMRG) method which has proved to be the best numerical tool for one-dimensional 
spin systems |15|. Ground state energy per site, the spin excitation gap and the two-spin 
correlation functions obtained from this method have been found to be accurate to several 
decimal places, in the case of the spin-| Heisenberg antiferromagnet which is amenable 
to exact study. In the DMRG method, spin parity symmetry can be used to characterize 
the spin states along with the s* ot as the good quantum numbers flEfl . There have been 
many interesting studies of spin-| and spin-1 chains in recent years, employing the DMRG 
technique [p^j . The spin-| and spin-1 dimerized Heisenberg chains with nearest neighbour 
and next nearest neighbour antiferromagnetic interactions have been recently studied by 
us by this new and powerful technique ||18||. There have also been some recent studies on a 
system with one or more spin-| impurities embedded randomly in a spin-1 chain[^], and 
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solution of spin models dynamically coupled to dispersionless phonons|2l| by the DMRG 
method. 

The spin wave analysis, in this paper, is followed by a report of our results from 
extensive DMRG studies on the alternating spin system. The DMRG calculations have 
been carried out on chains/rings with alternate spin-1 and spin-- sites. Studies on the 
ground state and low-lying excited states are reported in detail. Furthermore, by resorting 
to full diagonalization of the DMRG Hamiltonian matrix in different total M s sectors, we 
have also obtained the low-temperature thermodynamic properties of the alternating spin 
chain with periodic boundary condition. The thermodynamic properties we discuss will 
include low and high field magnetization, magnetic susceptibility and specific heat at low 
temperatures. These properties are compared with experimental studies on bimetallic 
chains. 

The paper is organized as follows. In the next section, we present properties of the 
ground and low-lying excited states obtained from spin wave analysis and DMRG calcu- 
lations on long alternating spin chain with and without dimerization. We also show that 
the short correlation length allows the ground state to be approximated by a variational 
ansatz. In the third section, we present the low-temperature thermodynamic properties 
of the alternating spin chain. We end the paper with a summary of all the results. 



2 Properties of the Ground and Low-Lying Excited 
states 

2.1 Spin- Wave analysis 

We begin with the Hamiltonian for a chain with spins s± and S2 on alternating sites. 

H = J y^[Sx, n • §2,71 + §2,71 • §1,71+1] (1) 

n 

where the total number of sites ( or bonds) is 2N, and we use periodic boundary conditions 
in this section with §i,iv+i = §i,i- 

The notation is illustrated in fig.l. We assume that si > S2, and will use spin wave 
theory [^TJ to compute the leading order corrections to the state shown in fig.l, where, 
the ^-component of the spin is si for the spin-si sites and — s 2 for spin-s 2 sites. The 
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Holstein-Primakoff transformations take the form 

l,n ~~ S l — °n a n 



for the spin-si sites, and 



= 2S 1 - aia n jan 

Kn = &i(^2 Sl -ala^j (2) 
S+ n = 6t^ 2 -^ 



S 2 " n = ^V2S2-fej6n (3) 

for the spin-S2 sites. We then rewrite the Hamiltonian in terms of the bosonic operators 
a n , b n , a|j and ¥ n , expand to quadratic order and Fourier transform to get 

H = -2NJ Sl s 2 + 2J^[s 1 6 t _ fe L fc + s 2 a\a k + ^/^cos(k/2)(a k b- k + d\b ] _ k )} (4) 

k 

This can be transformed using the Bogoliubov transformation, in the form 

Cfc = a k cosh 9 k + b < _ k sinh 
dfe = b- k cosh ^fc + a£ sinh 9 k 

tanh(2^,) = ^^cos(A;/2) (5) 
Si + s 2 

We then get 

if = -2NJ.Sx.S2 + Xl[^ifc4cfe + u 2k d\d k + (6) 
where the mode energies u;^ and zero-point energy z k are given by 

^lfc = J(-Si + s 2 )+w fc 

^fc = -J(si + s 2 )+^fc 



w fc = J v / ( s i- s 2) 2 + 4sis 2 sin 2 (fc/2) (7) 



5 



The ground state j^o > is the state annihilated by all the operators dk and dk, because 
ujik and U2k are positive for all k. We see that the modes denoted by ujik are gapless at 
k = 0, where they have a ferromagnetic dispersion uik ~ k 2 . The modes u 2 k are gapped 
for all k, with a minimum gap A = 2J(s\ — s 2 ) at k = 0. 

The ground state energy per bond is given by 

7T 

Eo t 1 [dk, T , x , , nN 

60 = 2iV = ~ JsiS2+ 2 J IP ' Sl + ^ + ^ O 



For si = 1 and s 2 = \, we get e = — 0.718J. The spin wave ground state \if) > can also 
be shown to be an eigenfunction of the S^ ot operator with the eigenvalue N(si — s 2 ). The 
sublattice magnetizations are given by the expectation values of S^ n and S^, 

1 1 f dk (si + s 2 )J 



o 2 / 1. 1 [dk 

<Sl, n > = (s 1 + -)--J- 







< S2,„ > = s 1 -s 2 -< S* n > (9) 
For si = 1, s 2 = L we find < >= 0.695 and < S z 2n >= -0.195. The operator 



S t + , = EK(V^) + fin(v^)] = c ( V2iV( Sl - s 2 ), (10) 

n 

to linear order in the bosonic operators. Since this annihilates \ip Q >, which is an eigenstate 
of Sf ot , we conclude that \ip > has s tot = s z ot = N(s 1 — s 2 ). Thus the ground state is a 
ferrimagnet. 

We can show that the u>ik modes are created by acting on \ipQ > with = -7= I](Sr,n+ 

ViV n ' 

S 2yn )e' ,kn , where k 7^ 0. The resultant states have s z ot = N(si — s 2 ) — 1 and are also 
annihilated by S^, t . Therefore we conclude that these states have Stot. = s lot. 
s 2 ) — 1. Since u>ik is a gapless branch, we further conclude that the system has gapless 
excitations to states with spin N(si — s 2 ) — 1. Similarly, the u 2 k modes are created by 
acting on \ip > with S^" = -7= XXS^ n + S 2n )e~ lkn , where k 7^ 0. By a similar argument, 

ViV n ' 

we can show that these states have s to t. — s z ot = N(si — s 2 ) + 1. The branch u 2 k is 
separated from the ground state by a gap A = 2J{s\ — s 2 ). Fig.2 shows these two 
excitation branches for s\ = 1 and s 2 — \ 

Finally, we can calculate the two-spin correlation functions. There are three kinds of 
functions one can consider, namely < Si i0 • Si >n >, < S 2 ,o • ^>\,n > an d < S 2 ,o • S 2 , n >. 
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We will consider only the first case as an illustrative example. Since the ground state has 
long range order ( with s to t. = s^ot. = N(s\ — s 2 ) ), we consider the subtracted correlation 
function 

< Si >0 • Si_ n > - < Si i0 > • < Si >n > = si < i) a \d a ] n + %a n \i/; > 

dk , , cos kn 



r dk , \ cos kn ,„ „ 

J 81 / — (si + s 2 ) (11 

J 7T OJh 




For n — > oo, we can show that the leading behaviour of this is given by an exponentially 
decaying factor e~ n ^, where £ _1 = m(si/s 2 ). For si = 1, s 2 — |, we get £ = 1.44. This 
remarkably short correlation length agrees well with numerical results as will be seen later. 
It may also be compared to the much larger values for the pure spin-1 antiferromagnet in 
which £ w 6|22], ^4| and the pure spin-| antiferromagnet in which £ = oo. 

We can use spin wave theory to study a dimerized model described by the Hamiltonian 
H = J£[(l + 8)S hn ■ S 2 , n + (1 - 5)S 2 , n • S ljn+1 ] (12) 

n 

where the dimerization parameter 5 lies in the range [0,1]. We find that the ground 
state and low-energy excitations are qualitatively similar to the undimerized case, 5 = 0. 
Namely, the ground state has Stot. = N(si — s 2 ). There is a gapless branch of excitations 
ui k with Stot. = N(si — s 2 ) — 1 and dispersion 



oJik = J(si + s 2 ) + Jy (si - s 2 ) 2 + 4sis 2 (l - 5 2 ) sm 2 (k/2), (13) 

and a gapped branch of excitations u 2 k with s tot . = N(si — s 2 ) + 1 and whose dispersion 
is, 

uJ 2 k = uJi k + 2 J (si - s 2 ) (14) 

To this order in spin wave theory, the gap at k = is given by A = 2J(si — s 2 ) independent 
of 5. Numerically, however we will see below that the gap increases almost linearly with 
5 for the spin-1 / spin-| alternating spin chain. 

2.2 Variational Calculation 

The short correlation length found above for the alternating spin-1 / spin-| chain suggests 
that a product wave function |22j could be a good variational trial state for obtaining the 
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ground state properties. We work in the Fock space basis of the alternating spin chain. We 
are interested in a variational wave function in which the state |... 1 — | 1 —\ . . . > 
of fig. 1 is expected to have the largest amplitude, say, 1. We can then show that a state 
in which there are n spins on the spin-1 sublattice in the state s z = must have a real 
amplitude with the sign (—1)™. This sign rule (which is analogous to the Marshall sign 
rule for the spin-| chain J23|) can be proved by using the Perron- Frobenius theorem and 
the fact that the states with an odd number of spins in the s z = state are connected 
to states with an even number of spins in the s z = states by the operation of an odd 
number of exchange terms in the Hamiltonian (Eq.l). 

Let us now introduce the following real and positive amplitudes for the six states 
possible for a nearest neighbour bond, namely, unit amplitude for |1, — ~ > and | — 1, | >, 
an amplitude rji for |0, — | > and |0, \ > and amplitude r/ 2 for | — 1, — | > and |1, | >. 
Now, we consider a variational wave function of the form 

l*fai,%)>=£C a |V«> (15) 

a 

where a runs over all the 6^ possible states of N spin-1 and spin-^ sites and C a = 

(— l) n T7^a,j, where the product is over all the 2N bonds, u a j is the amplitude of the 

j 

bond j and n denotes the number of spins in the s z = state in the chain state a. 

The state |\l/(77i, 772) > is translationally invariant. However, it is not an eigenstate 
of the total spin Sf ot or even Sf ot ; further, since it includes states {s^} and {— s z n } with 
equal amplitude, the expectation value of Sf n or Sf n is zero for any site n. In spite of 
these drawbacks, we will see below that the 1^/(771,772) > gives a good variational ground 
state energy for an appropriate choice of 7/1 and 772. 

Calculations involving the state 1^/(771,772) > can be carried out using the transfer 
matrix method. For instance, the probabilities of the six possible bonds for \\,n > < 
l,n + 1| are given by the 2x3 matrix. 

A = ( V2 vl 1 

v 1 v! vi 

while the probabilities of the bonds |1, n > (g) < |, n\ are given by the 3x2 matrix A T . 
Then the norm < ^(771, 77 2 )|\I/(r/i, 772) >= Tr(AA T ) N . The matrix AA T has the eigenvalues 
A + = (1 + ^)2,^2^4 and ^_ _ ^_^2)2 Hence, < ^(771,772)1^(771,772) >= A^ + A^, which 
is dominated by as N — > 00. 
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We can now calculate the variational energy per bond e(r]i, 772) through the expression 

e(^,%) < tf(77i,772)|¥(?7i,»k) >= ^ < ^2) 1^1*^1,772) > (16) 

This is equal to < ^(r]i,r] 2 )\Si n ■ Si >ri+ i |\I/(?7i, 772) > by translational invariance. We find 
that, 

[-1(1 + r/ 2 2 ) 3 (l - rj) - 2yg(l + 7^)^(1 + V2 ) - ^(1 - gj) - 2V2qg(l + 77 2 )] 

[(l + 7ji) 2 + 2^p 1 <} 

We now minimize this function of r/i and 772, and find that the minimum occurs at 771 = 
0.842, 772 = 0.445, giving e = — 0.701J. This compares favourably with the spin wave 
result of -0.718J. 

We can compute the two-spin correlation function and determine how it decays asym- 
potically at large distances. The correlation length £ for the asymptotic decay is given 
by, £~ 1 = m(A_|_/A_). Hence £ is 0.749 at the values of 771, 772 given above. This is even 
shorter than the value of 1.44 found from the spin wave theory. 

We can improve the variational calculation by allowing the five independent ( taking 
the amplitude for ~ >= 1) real amplitudes for the different possible states of the 
nearest neighbour bonds rather than only two amplitudes ( 771, t? 2 ) as used above. However 
we will not do so here as such a calculation would be tedious and the two-amplitude ansatz 
has already given us a good understanding of the short correlation length and the ground 
state energy. 



2.3 DMRG Studies 

We have performed DMRG calculations on alternating spin-1 / spin-| chain with open 
boundary condition for the Hamiltonian in Eq.(l). We compute the ground state proper- 
ties by studying chains of upto 50 to 100 sites. The number of dominant density matrix 
eigenstates, m, that we have retained at each DMRG iteration is between 80 to 100. The 



DMRG procedure follows the usual steps for chains discussed in earlier papers ]T5|, |T7[ [I8| 
except that the chains do not have the symmetry between the left and right halves and 
the density matrices for these two halves have to be constructed at every iteration of the 
calculations. The DMRG results reported in this paper are all specialized to the case of 
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si = 1 and S2 = \- The ground state of the chain lies in the M s = N(si — S2) sector for 
a 2N sites system and we cannot use spin parity symmetry in this sector. In general, for 
M s 7^ 0, parity cannot be exploited, keeping M s a good quantum number. The ground 
state is confirmed to be in M s = N(s\ — S2) sector for a chain of 2N sites from extensive 
checks carried out by obtaining the low-energy eigenstates in different M s sectors of a 
20-site chain. A state corresponding to the lowest energy in M s = N(s\ — S2) is found in 
all subspaces with \M S \ < N{s\ — $2) and is absent in subspaces with M s > N(si — S2). 
This shows that the spin in the ground state, sg = N(s\ — S2). In fig. 3, we show the 
extrapolation of the energy per site as function of inverse system size. The energy/site, 
eo, extrapolates to — 0.72704J for the ground state of the system. The spin wave analysis 
value of — 0.718J for eo is higher than the DMRG value. It is worth mentioning at this 
point that the DMRG ground state energy/site agrees to better than 1CT 7 with the exact 
solution in the case of uniform spin-| Heisenberg antiferromagnet JT5| It is interesting 
to note that, in the alternating spin case, the energy/site lies in between the values for 
the pure spin-| uniform chain ( — 0.443147J ) and the pure spin-1 uniform chain ( 
-1.401484J ) ||. 

In fig. 4, we show the spin-density at all the sites of a chain of 100 sites. The spin- 
density is uniform on each of the sublattices in the chain. The expectation value of S z at 
the spin-1 site reduces from the classical value of 1 to 0.79248, while, at the spin-^ site, it is 
—0.29248. This can be compared with spin wave values of 0.695 and —0.195 for spin-1 and 
spin- \ sites respectively. We note that the spin wave analysis overestimates the quantum 
fluctuations. It is very interesting to note that, the spin-density distribution in alternating 
spin chain behaves more like in a ferromagnetic chain than like in an antiferromagnet, 
with the net spin of each unit cell perfectly aligned ( but with small fluctuations on 
the individual sublattices). In a ferromagnetic ground state the spin-density at each site 
has the classical value appropriate to the site spin, whereas for an antiferromagnet, this 
averages out to zero at each site as the ground state is nonmagnetic. From this viewpoint, 
the ferrimagnet is similar to a ferromagnet and is quite unlike an antiferromagnet. The 
spin wave analysis also yields the same physical picture. 

Owing to the alternation of the spin-1 and spin-| sites along the chain, one has to 
distinguish between < S^ S^ n >, < S^o^n > an d < ^io^2n > P & i r correlations. We 
calculate all the three correlation functions with the mean values subtracted out as in 
Eq.(ll), since the mean values are nonzero in this system unlike in a pure antiferromag- 
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netic spin chain. In the DMRG procedure, we have computed these correlation functions 
from the sites inserted at the last iteration, to minimize numerical errors. In fig. 5, we 
plot the two-spin correlation functions in the ground state as a function of the distance 
between the spins for an open chain of 100 sites. All three correlation functions decay 
rapidly with distance. From the figure it is clear that, except for < Sf S| n > corre- 
lation, other correlations are almost zero even for the shortest possible distances. The 
< S* Q S| n > correlation has a finite value ( —0.094 ) only for the nearest neighbours. 
This rapid decay of the correlation functions do not easily allow one to find the exact cor- 
relation length £ for a lattice model though it is clear that £ is very small (i.e. less than 
one unit). As mentioned earlier, spin wave theory gives £ = 1.44 while the variational 
calculation gives £ = 0.75. This type of decay could be compared with the behaviour of 
the correlation function at the Majumdar- Ghosh point [^6) for the pure spin-| Heisenberg 
chain with nearest and next nearest neighbour antiferromagnetic exchange interactions 
( J\ — J 2 model at J 2 = 0.5Ji). The AKLT model at the exactly solvable point 



in 



the case of a pure spin-1 chain described by a bilinear-biquadratic Hamiltonian also has a 
very short correlation length ( £ = 0.91 ). Both the cases compared above, however, have 
one property in common, i.e., the exactly solvable point in both models lie in the gapped 
phase. 

The lowest spin excitation of the alternating spin chain is to a state with s = sq — 1. 
To get this state, we target the 2nd state in M s = sg — 1 sector of the chain. To confirm 
that this state is the s = sq — 1 state, we have computed the 2nd state in M s = sector 
and find that it also has the same energy. However, the corresponding state is absent in 
M s sectors with M s > \s@ — 1|. Besides, from exact diagonalization of all the states of an 
alternating spin chain with 8 sites, we find that the ordering of the state is such that the 
lowest excitation is to the state with spin s = sg — 1- We have obtained the excitation gap 
of the alternating spin chain in the limit of infinite chain length by extrapolating from the 
plot of spin gap vs the inverse of the chain length ( fig.6 ). We find that this excitation 
is gap less in the infinite chain limit. This is very unusual in that the correlation function 
in the ground state decays exponentially ( in fact £ < 1) but the system has a gap less 
excitation, in agreement with the spin wave analysis. 

Motivated by the spin-wave calculation, we also have computed the energy of the 
s = sg + 1 state by targetting the lowest state in M s = sg + 1 sector. In fig. 7, we plot 
the excitation gap to the s = sg + 1 state from the ground state as a function of the 
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inverse of the chain length. The gap saturates to a finite value of (1.2795 ± 0.0001) J. 
The S z expectation values computed in this state are found to be uniform on each of the 
sublattices ( independent of the site) and leading us to believe that this excitation is not 
a magnon like excitation ( quantized in a box ) as has been observed for a spin-1 chain 
in the Haldane phase fTCf] . 

It would be interesting to know the total spin of the states as a function of their 
energies. For a smaller alternating spin system, it is possible to characterize all the states 
by their energy and total spin value, by resorting to an exact diagonalization scheme. 
The total spin value of a state is naturally fixed if we exploit the total spin conservation 
property of the Hamiltonian while constructing the Hamiltonian matrix. This can, for 
example, achieved by using a valence bond basis|2^] for setting up the Hamiltonian matrix. 
Alternately, we can also compute the expectation value of the total spin operator in each 
eigenstate in the M s = sector to provide a spin label for each of the states. We have 
followed the latter procedure. In fig. 8, we present the energy levels as a function of 
the total spin of the states for an eight site ferrimagnetic ring and twelve site spin-^ 
ferromagnetic and antiferromagnetic rings. We find that the spin of the state appears to 
vary irregularly with energy unlike in the case of pure spin-- ferro and antiferromagnets 
(fig. 8). Careful comparison of the three plots in the figure shows that the low-lying 
excitations of the ferrimagnet to spin states s to t < $g is ferromagnetic like and to states 
with Stot > s g is antiferromagnetic like, for finite systems. 

We have also studied the dimerized alternating spin chain defined by the Hamiltonian 
( Eq.12 ) with 5, the dimerization parameter, in the range < 8 < 1. Earlier works 



on spin chains |28|] have revealed that with the alternation 8 in the exchange parameter, 
the half integer spin chain will have an unconditional spin-Peierls transition whereas for 
integer spin chain the transition is conditional. This conclusion has been drawn from the 
fact that, with the inclusion of 5, the magnetic energy gain AE can be defined as, 

AE(2N, 8) = ^[E(2N, 5) - E(2N, 0)] (18) 

where E(2N, 5) is the ground state energy of the 2N sites system with alternation 5 in 
the exchange integral and E(2N, 0) is the ground state energy of the uniform chain of 
2N sites. For the pure spin chain, if we assume that AE varies as 5 U for small 5, we find 
that v = 4/3 for the spin-| chain and v = 2 for the spin-1 chain p8|. Thus, for the spin-| 
chain, the stabilization energy always overcomes the elastic energy, whereas for the spin-1 
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case, it depends on the stiffness of the lattice. 

We have employed a DMRG calculations to obtain AE(2N,S), for small values of S 
for the alternating spin chain. To determine the exact functional form of the magnetic 
energy gain, we varied the chain length from 50 sites to 100 sites and also m( the number 
of states retained in each DMRG iteration) from 80 to 100 to check the convergence of 
AE with the chain length. The dependence of AE(2N,5) on 1/2N is linear for the 5 
values we have studied. Fig. 9 gives a sample variation of AE(2N,5) upon 1/2N. This 
allows us to extrapolate AE(2N,5) to the infinite chain limit reliably. In fig. 10, we show 
the plot of AE(2N, 5) vs. 5 for finite 2iV values and also the extrapolated infinite chain. 
We see that there is a gain in magnetic energy upon dimerization even in the infinite 
chain limit. To obtain the exponent is, we plot In AE(2N, 5) vs. ln<5 for the infinite chain 
( fig. 11 ). From this figure, we find that in the alternating spin case, for the infinite chain 
AE ~ S 2m±om . Thus, the spin-Peierls transition appears to be close to being conditional 
in this system. The magnetic energy gain/site for finite chains is larger than that of the 
infinite chain for any 5 values ( fig. 10). It is possible that the distortion in finite chain is 
unconditional while that of the infinite chain is conditional. 

We have also studied the spin excitations in the dimerized alternating spin-1 / spin-| 
chain. We calculate the lowest spin excitation to the s = s G — 1 state from the ground 
state. We find that the s = sq — 1 state is gapless from the ground state for all 5 values. 
This result agrees with the spin wave analysis. The system remains gapless even while 
dimerized unlike the pure antiferromagnetic dimerized spin chains. There is a smooth 
increase of the spin excitation gap to s = sq + 1 state from ground state with increasing 
5. We have plotted this gap with 5 in fig. 12. The gap shows almost a linear behaviour as 
a function of 5, with an exponent of 1.07 ± 0.01. The spin wave analysis however shows 
that this excitation gap is independent of 5. 

3 Low- Temperature Properties 

In this section, we present results of our DMRG calculations of the thermodynamic prop- 
erties of the alternating spin system. The size of the system varies from 8 to 20 sites. 
We impose periodic boundary conditions to minimize finite size effects. We set-up the 
Hamiltonian matrices in the DMRG basis for all allowed M s sectors for a ring of 2N sites. 
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We can diagonalize these matrices completely to obtain all the eigenvalues in each of the 
M s sectors. As the number of DMRG basis states increases rapidly with increasing m, we 
retain a smaller number of dominant density matrix eigenvectors in the DMRG procedure, 
i.e., 50 < m < 65, depending on the M s sector as well as the size of the system. We have 
checked to find the dependence of properties ( with m in the range 50 < m < 65 ) for 
the system sizes we have studied ( 8 < 2N < 20 ) and have confirmed that the properties 
do not vary significantly for the temperatures at which they are computed. The above 
extension of the DMRG procedure is found to be accurate by comparing with exact diag- 
onalization results for small systems. It may appear surprising that the DMRG technique 
which essentially targets a single state, usually the lowest energy state in a chosen sector, 
should provide accurate thermodynamic properties since these properties are governed by 
energy level spacings and not the absolute energy of the ground state. However, there 
are two reasons why the DMRG procedure yields reasonable thermodynamic properties. 
Firstly, the projection of the low-lying excited state eigenfunctions on the DMRG space 
in which the ground state is obtained is substantial and hence these excited states are well 
described in the chosen DMRG space. The second reason is that the low-lying excitations 
of the full system are often lowest energy states in different sectors in DMRG procedure 
and thus their energies are quite accurate even on an absolute scale. 

The canonical partition function Z for the 2N site ring can be written as 

z = ^ e -m-B(M s)j ) (19) 
j 

where, the sum is over all the DMRG energy levels of the 2N site system in all the M s 
sectors. Ej and (M s )j are energy and z-component of the total spin of the state j, B is the 
strength of the magnetic field in units of J/gfiB ( 9 is the gyro magnetic ratio and \ib is the 
Bohr magneton ) along z direction and (3 = J/k B T with k B and T being the Boltzmann 
constant and temperature respectively. The field induced magnetization, < M >, can be 
defined as 

T(M) .e-PiEj-BiM^) 
< M >= ^ Ms) i e (20) 

Zj 

the magnetic susceptibility, Xi by relating it to the fluctuation in magnetization, 

X = P[< M 2 > - < M > 2 ] (21) 
and similarly the specific heat, C, by relating it to the fluctuation in energy, can be written 
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as, 

C = ^[<E 2 >-<E> 2 ] (22) 

The dimensionalities of the DMRG Hamiltonian matrices that we completely diag- 
onalize vary from 2500 to 3000, depending upon the DMRG parameter m and the M s 
value of the targetted sector, for rings of sizes greater than 12. These matrices are not 
very sparse, owing to the cyclic boundary condition imposed on the system. The DMRG 
properties compare very well with exact results for small system sizes amenable to exact 
diagonalization studies. In the discussion to follow, we present results on the 20-site ring 
although all calculations have been carried out for system sizes from 8 to 20 sites. This 
is because the qualitative behaviour of the properties we have studied are similar for all 
the ring sizes in this range. 

We present the dependence of magnetization on temperature for different magnetic 
field strengths in fig. 13. At low magnetic fields, the magnetization shows a sharp decrease 
at low-temperatures and shows paramagnetic behaviour at high temperatures. As the field 
strength is increased, the magnetization shows a slower decrease with temperature and 
for field strength comparable to the exchange constant, the magnetization shows a broad 
maximum. This behaviour could be understood from the type of spin excitations present 
in the system. The lowest energy excitation at low magnetic fields is to a state with spin 
s less than so- Therefore, the magnetization initially decreases at low temperatures. As 
the field strength is increased, the gap to spin states with s > sq decreases as the Zeeman 
coupling to these states is stronger than to the states with s < sq- The behaviour of 
the system at even stronger fields turns out to be remarkable. The magnetization in the 
ground state ( T = ) shows an abrupt increase signalling that the ground state at this 
field strength has M s > s G . The temperature dependence of the magnetization shows a 
broad maximum indicating the presence of states with even higher spin values lying above 
the ground state in the presence of this strong field. Only at very intense fields do we find 
the magnetization decreasing with increasing temperature. This happens because at such 
large field strengths, the ground state is the highest spin state possible for the system. 

The dependence of magnetization on the magnetic field is shown at different temper- 
atures in fig. 14. At low temperature the magnetization shows a plateau. The width of 
the plateau decreases as the temperature is raised and eventually the plateau disappears. 
The existence of the plateau shows that the higher spin states are not accessible at the 
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chosen temperature. At higher fields, the larger Zeeman splittings of higher spin states 
become accessible leading to an increase in the magnetization. All these curves intersect 
at B — J I qhb and B = 2.5J/gp,B and these fields are close to the field strengths at which 
the ground state switches from one M s value to another higher value. 

The dependence of xl '/2N on temperature for different field strengths are shown in 
fig. 15. For zero field, the zero temperature value of xT is infinite in the thermodynamic 
limit and for finite rings is finite and equal to the fluctuation in magnetization. For the 
ferrimagnetic ground state xT/2N, as T — > 0, is given by sq(sg + l)/6N. As the tem- 
perature increases, this product decreases and shows a minimum around k^T = 0.5J 
before increasing again. The minimum manifests due to the states with M s < sq getting 
populated at low-temperatures. These states in the infinite chain limit turn out to be 
the gapless excitations in the spectrum. The subsequent increase in the %T product is 
due to the higher energy and higher spin states accessed with further increase in tem- 
perature. Experimentally, it has been found in the bimetallic chain compounds that the 
temperature at which the minimum occurs in the \T product depends upon the magni- 
tude of the spins s\ and |P-3| . The Ni 11 — Cu 11 bimetallic chain shows a minimum in 
XT/2N at a temperature corresponding to 55 cvrC 1 (80 K) and independent estimate of 
the exchange constant in this system is 100 cm' 1 129|]. This is in very good agreement 
with the minimum theoritically found at temperature (0.5 ± 0.1)J. The minimum in 
XT/2N vanishes at B — O.U/gfiB corresponding to ~ 10T and it would be interesting to 
study the magnetic susceptibility of other systems experimentally under such high fields. 
The low-temperature zero-field behaviour of our system can be compared with the one- 
dimensional ferromagnet. In the latter, the spin wave analysis shows that the \T product 
increases as ^ at low temperatures pOfl . 

In finite but weak field, the behaviour of xT is different. The magnetic field opens up 
a gap and xT goes exponentially to zero for temperatures less than the gap in the applied 
field. Even in this case a minimum is found at the same temperature as in the zero-field 
case, for the same reason discussed in the zero field case. 

In stronger magnetic fields, the behaviour of xT from zero temperature upto k^T — 
0.5 J is qualitatively different. The minimum in this case vanishes. In these field strengths, 
the states with higher M s values are accessed even below ksT = 0.5 J. The dependence 
of xT above k^T = 0.5J is the same in all cases. In even stronger magnetic fields, the 
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initial sharp increase is suppressed. At very low temperature, the yT product is nearly 
zero and increases linearly with T over the temperature range we have studied. This 
can be attributed to a switch in the ground state at this field strength. The very high 
temperature behaviour of yT should be independent of field strength and should saturate 
to the Curie law value corresponding to the mean of spin-1 and spin-| values which is 
11/24. 

The temperature dependence of specific heat also shows marked dependence on the 
magnetic field at strong-fields. This dependence is shown in fig. 16 for various field 
strengths. In zero and weak magnetic fields, the specific heat shows a broad maximum 
around /c#T = 0.6J. At strong-magnetic field (B = J), there is a dramamtic increase in 
the peak height at about the same temperature, although the qualitative dependence is 
still the same as at low-magnetic fields. This indicates that the higher energy high-spin 
states are brought to within ksT of the ground state at this magnetic field strength. 

Studies on dimerized alternating spin chains show qualitatively similar trends as the 
uniform chains. This is not surprising as the low-energy spectrum of the system does not 
change qualitatively upon dimerization. 

4 Summary 

We have studied the alternating spin-1 / spin-| model in detail. The ground and low- 
lying excited states have been analyzed by using a spin wave theory as well as DMRG 
calculations. Both the methods predict a ground state with spin s G = N(s 1 —s 2 ) for a 2N 
site system. They also predict a gapless excitation to a state with s = sq — 1 in the infinite 
chain limit. The lowest gapped excitations are to states with spin s — sq + 1- The very 
short correlation length in the ground state of the system motivated its description by a 
variational trial function of the product type. Interestingly, the spectrum is qualitatively 
unchanged upon dimerization. The dimerization is itself conditional in the infinite chain 
limit. 

The DMRG technique is also employed to obtain the low-temperature thermodynamic 
properties. The magnetic susceptibility shows very interesting magnetic field dependence. 
The yT vs. T plot shows a minimum at low-magnetic fields. This minimum vanishes at 
high-magnetic fields. The specific heat shows a maximum as a function of temperature 
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at all fields. The height of the maximum shows a dramatic increase at high-magnetic 
field. Experimental systems describable by this model exist and have been studied quite 
extensively. It is hoped that our studies will motivate experimental studies of these 
systems in high magnetic fields. 
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Figure Captions 

Fig.l 

Schematic picture of the arrangement of spins s\ and s 2 on a chain with interactions 
discussed by the Hamiltonian in Eq.(l). 

Fig.2 

The two branches of the spin wave dispersion curves, uj^ and u 2 k ( Eq.(7) ), in units of 
J, for the alternating spin system with si = 1 and s 2 — \- 

Fig.3 

Extrapolation of the ground state energy/site (eo), in units of J, as a function of inverse 
system size, for the uniform alternating spin-1 /spin-| chain. 

Fig.4 

Expectation value of the z-component of the site spin vs. site index, n. The upper and 
lower points are for the spin-1 (i = 1) and the spin-| (i = 2) sites respectively. 

Fig.5 

Subtracted two-spin correlation functions (defined in the text) as a function of distance 
between the two spins, (a) spin-1 spin-1 correlations, (b) spin-| spin-| correlations and 
(c) spin-1 spin-| correlations. 

Fig.6 

Energy gap (units of J) from the ground state to the lowest state with spin s = s G — 1 
as a function of inverse system size, sg is the total spin of the ground state. 

Fig.7 

Excitation gap (units of J) from the ground state (spin s = sg) to the state with spin 
s = sg + 1 as a function of the inverse system size. The convergence to the infinite system 
is much faster for this gapped excitation, as compared to the gapless excitation in fig.6. 

Fig.8 

Plot of energy (units of J) vs. the total spin quantum number of the complete spectrum 
of (a) an 8 site ring with an alternating spin-1 / spin-| arrangement (b) a ring of 12 site 
spin-i antiferromagnet and (c) a ring of 12 site spin-| ferromagnet. The states with very 
high energies are not shown, for (b) and (c). 

Fig.9 

Gain in magnetic energy (units of J) associated with dimerization vs. the inverse system 
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size for three different values of dimerization 8. (a) 5 = 0.025, (b) 5 = 0.05 and (c) 
5 = 0.075. 

Fig.10 

Magnetic energy gain AE(2N,5) (units of J) as a function of dimerization parameter 
S for different system sizes. In the figure 2N = 50 (squares), 2N = 100 (circles) and 
extrapolated values with N — > oo (triangles) are shown. 

Fig. 11 

Log-log plot of extrapolated magnetic energy gain (units of J) for infinite system size and 
dimerization parameter 5. The slope is calculated to be 2.00 ± 0.01. 

Fig.12 

Excitation gap (units of J) to the state with spin s = sq + 1 from the ground state 
(s = sq) as a function of S for the dimerized alternating spin-1 / spin-| chain. 

Fig.13 

Plot of magnetization/site as a function of temperature T for five different values of 
magnetic fields, B. Squares are for B = 0.1J/gfi B , circles for B = 0.5J/gfi B , triangles for 
B = J/g[i B , diamonds for B = 2J/gfi B and inverse triangles for B = ?>J/gn B . 

Fig.14 

Magnetization/site vs. the magnetic field strength B, in units of J/gfi B , for four different 
temperatures T. T = 0.3 J/ k B results are given by squares, T = 0.5 J/ k B by circles, 
T = 0.7 J I k B by triangles and T = J/k B by diamonds. 

Fig.15 

XT ( defined in the text) per site as a function of temperature T for various magnetic 
field strengths, B. Zero field results are shown by squares, B = 0.01J/gfi B by circles, 
B = 0.1J/gfi B by triangles and B = J j g\i B by diamonds. 

Fig.16 

Specific heat/site as a function of temperature T for four different values of magnetic 
fields, B. Zero field data are shown by squares, B = 0.01J/gfi B by circles, B = 0.1J/gfj, B 
by triangles and B = J jg\i B by diamonds. 
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